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Abstract 

We develop efficient algorithms for level-truncation computations in open bosonic string field 
theory. We determine the classical action in the universal subspace to level (18,54) and apply 
this knowledge to numerical evaluations of the tachyon condensate string field. We obtain two 
main sets of results. First, we directly compute the solutions up to level L = 18 by extremizing 
the level-truncated action. Second, we obtain predictions for the solutions for L > 18 from an 
extrapolation to higher levels of the functional form of the tachyon effective action. We find that 
the energy of the stable vacuum overshoots -1 (in units of the brane tension) at L = 14, reaches 
a minimum Emin = —1.00063 at L ~ 28 and approaches with spectacular accuracy the predicted 
answer of -1 as L —> oo. Our data are entirely consistent with the recent perturbative analysis 
of Taylor and strongly support the idea that level-truncation is a convergent approximation 
scheme. We also check systematically that our numerical solution, which obeys the Siegel gauge 
condition, actually satisfies the full gauge-invariant equations of motion. Finally we investigate 
the presence of analytic patterns in the coefficients of the tachyon string field, which we are able 
to reliably estimate in the L —> oo limit. 
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1 Introduction and Summary 


The realization that D-branes are solitons of the open string tachyon [||, ^ has triggered a revival 
of interest in open string field theory. Much work has focused on the search for classical solutions 
of cubic bosonic open string field theory Q (OSFT). Despite important technical progress in the 
understanding of the open string star product - notably the discovery of star algebra projectors 
1^, ^ and of new connections with non-commutative field theory I, 0, I, i - analytic classical 
solutions of OSFT are still missing. 


Fortunately, the OSFT equations of motion can be solved numerically in the ‘level-truncation’ 
scheme invented by Kostelecky and Samuel |jT^. The open string field is restricted to modes 
with an Lq eigenvalue smaller than a prescribed maximum ‘level’ L. For any hnite L, the trun¬ 
cated OSFT action contains a finite number of fields and numerical computations are possible. 
Remarkably, numerical results |17|-|37] for various classical solutions appear to converge rapidly 
to the expected answers as the level L is increased. Much of our present intuition about the 
classical dynamics of OSFT comes from the level truncation scheme, and in fact even in vacuum 
string field theory 00 several exact results were first guessed based on numerical data. 


This motivated us to develop efficient algorithms for level-truncation calculations. Our main 
technical innovations are the systematic use of conservation laws to compute the cubic ver¬ 
tices, and the implementation of our algorithms on a C-|--|- code. In this paper we apply these 
methods to the evaluation of the classical action in the universal subspace P, which is the 
space of string fields generated by ghost oscillators and matter Virasoro generators acting on 
the vacuum. Using conservation laws we determine the classical action directly in the universal 
basis in a recursive way, with an algorithm whose complexity is linear in the number of vertices 
(cubic in the number of fields). Some details about the numerical algorithms can be found in 
appendix A of the paper. The gain in efficiency of our methods is of several orders of magnitude, 
and we are able to obtain the 10^*^ universal cubic vertices at level (18,54). 


The universal subspace has special physical significance because it contains the tachyon 


condensate string held, the solution of OSFT corresponding to the stable vacuum of the open 


string tachyon. Its (negative) energy per unit volume must exactly cancel the D-brane tension. 
Sen and Zwiebach’s computation |^| of the tachyon condensate up to level (4,8) gave the first 
evidence that OSFT reproduces the correct D-brane physics. Moeller and Taylor ||20|| pushed 
the computation to level (10,20) finding that 99.91% of the D-brane tension is cancelled in 


^ Exact results have been obtained in vacuum string field theory ||^, (VSFT), see e.g. 10 li-00- 
VSFT is the version of open string field theory which appears to describe the stable tachyon vacuum. In VSFT 
the BRST operator is replaced by a purely ghost insertion at the string midpoint and the equations of motion 
take the exactly solvable form of projector equations (with an auxiliary ‘twisted’ ghost system [hill. It would 
nevertheless be extremely desirable to solve analytically the original OSFT equations. See e.g. |p_q| for some 
formal attempts. 
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L 

E(L,3L) 

E{L,2L) 

2 

-0.9593766 

-0.9485534 

4 

-0.9878218 

-0.9864034 

6 

-0.9951771 

-0.9947727 

8 

-0.9979302 

-0.9977795 

10 

-0.9991825 

-0.9991161 

12 

-0.9998223 

-0.9997907 

14 

-1.0001737 

-1.0001580 

16 

-1.0003754 

-1.0003678 

18 

-1.0004937 

-1.00049 


Table 1: Values of the vacuum energy in level-truncation, in the {L,3L) and (L, 2L) approxi¬ 
mation schemes. 


the tachyon vacuum. Given such a remarkable agreement, it may appear quite pointless to 
extend their results to higher level. Not so. Up to level 10, the individual coefficients of the 
string field appear to converge much less rapidly than the value of the action. A more precise 
determination of the coefficients is likely to provide clues for an exact solution. Indeed various 
surprising patterns obeyed by OSFT solutions were ‘experimentally’ observed in [^|, but more 
accurate results are needed to decide which of these patterns are likely to be exact and which 
ones only approximate. Higher level computations can also be expected to shed light on the 
nature of the level truncation procedure itself, which still lacks a sound theoretical justification. 


Our first main set of results (described in section 3) is the computation of the Siegel gauge 
tachyon condensate in level-truncation up to L = 18. The procedure is the standard one: at 
any given level L, there are Nl scalar helds that obey the Siegel condition, and we determine 
their vev’s by solving the equations of motion implied by the gauge-fixed action. There 
is a potential subtlety here: the full equations of motion before gauge-fixing impose a bigger 
number of constraints [^| (the extra conditions simply enforce extremality of the action along 
gauge orbits). Consistency demands that the full set of equations of motion is satisfied as 
L —> oo, and we systematically check that this indeed happens, with remarkable accuracy. As 
another consistency check, we verify that the tachyon condensate obeys the quadratic relations 
analytically derived by Schnabl and we again hnd excellent agreement. 


The values of the vacuum energy as a function of L are shown in Table |I|. Unexpectedly, 
at L = 14 the energy overshoots the predicted answer of -1 and appears to further decrease 
at higher levels. As a first reaction, one may wonder whether the level-truncation procedure is 
breaking down for L > 10, as could happen if the approximation was only asymptotic. In this 
pessimistic scenario, for any OSFT observable there would be a maximum level that gives the 
estimate closest to the ‘exact’ value, and beyond this optimal level the procedure would stop 
converging. However, the data favor a smooth behavior as L increases, since the differences 
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between consecutive approximations are getting smaller. 

The results in Table 1 may simply indicate that the approach of the energy to -1 as L —> oo is 
non-monotonic, contrary to previous naive expectations. Indeed, Taylor has recently presented 
convincing evidence for this benign interpretation of our results.^ He applies a clever 
extrapolation technique to level-truncation data for L < 10 to estimate the vacuum energies 
even for L > 10. This procedure reproduces quite accurately our exact values in Table 1 and 
further predicts that the vacuum energy reaches a minimum for L ~ 28, but then turns back to 
approach asymptotically -1 for L ^ oo. 


E 



Figure 1: Curves of the vacuum energy as a function of level, as predicted by our extrapolation 
scheme for various values of M (maximum level of the data used in the extrapolation). The 
figure shows the curves E^^\L) on a logarithmic plot, for M between 8 (lowermost curve) and 
16 (uppermost curve). Data in the (L,3L) scheme. 

In sections 4 we introduce our second main set of results. We devise an extrapolation 
technique in the same spirit of Taylor’s analysis. We consider the effective tachyon potential 
Vl{T) around the unstable vacuum, obtained by classically integrating out all the higher scalars 
up to level L. Vl{T) is computed ‘non-perturbatively’ by fixing the value of T and solving 
^The data in Table 1 were first announced at the Strings 2002 conference, Cambridge, July 2002 0 . 
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Lmin 

E 

^min 

Eoo 

M = 6 

41.1 

-1.001171 

-1.000949 

M = 8 

28.2 

-1.000660 

-1.000140 

M = 10 

27.8 

-1.000646 

-1.000113 

M = 12 

27.5 

-1.000637 

-1.000077 

M = 14 

27.3 

-1.000633 

-1.000046 

M = 16 

27.3 

-1.000632 

-1.000030 


Table 2: Parameters of the curves (in the {L,3L) scheme). The energy reaches its 

minimum Emin for L = Lmin, and tends asymptotically to E^o as L —> oo. 

numerically the equations of motion for the other scalars^. We are able to obtain Vl{T) up to 
L = 16. Clearly for each L, the minimum of Vl{T) is just the vacuum energy E^ at level L. 
However the functional dependence on T contains more information than just the extremal value 
El- The idea is to perform an extrapolation in L of the whole functions in T. In practice, we 
consider a finite interval of values of T around the non-perturbative minimum. For a fixed T in 
this interval we interpolate our data for Vl{T) with a polynomial in 1/L, and then extrapolate 
this polynomial to higher levels. To check the stability of this approximation scheme, we vary 
the maximum level M of the set of data used as input for the extrapolation: for each M < 16, 
we apply the extrapolation to the functions {Vl{T) \ L < M}. This gives estimates for the 
tachyon vev and for the corresponding vacuum energy E^"^ , for any L > M. 

The predicted power of the method is quite impressive. For example, with M = 10, that is 
using only level-truncation results up to level 10, the estimate reproduces with an accuracy 
of 10“® the exact tachyon vev Tig, obtained by straightforward level-truncation at L = 18. This 
is remarkable, since the former computation is over a thousand times faster than the latter. 

Figure || and Table § summarize the extrapolations of the energy as a function of level, for 
various values of M. The data completely confirm (with enhanced precision) the conclusions of 
Taylor The behavior of the energy as a function of level is non-monotonic, but eventually 
the asymptotic limit of -1 is reached with spectacular accuracy. 

These results greatly reassure us of the validity of the level-truncation scheme. Observables 
in OSFT have a smooth limit as L —> oo, which (in the absence of an alternative definition) 
should be identified with their ‘exact’ value. In all cases where an independent prediction for 
the observable is available (as for the vacuum energy, or for Schnabl’s quadratic identities), the 
L —> oo extrapolation gives the correct answer. 

A practical lesson of this analysis is that polynomial interpolations in 1/L have great pre¬ 
dictive power, at least for the (L,3L) approximation scheme^. This observation makes the 

®In this we differ from Taylor who uses instead a series expansion of the potential in powers of T. 

In the (L, 3L) scheme, the string field is truncated up to level L and all of its mode are kept in the OSFT 
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level-truncation scheme much more efficient, as reliable estimates can be extracted from (rela¬ 
tively) painless numerical work. 

In section 5 we describe the results for the individual coefficients of the tachyon string field 
extrapolated to L = oo. The main conclusions are: 


• The asymptotic value of the tachyon vev is Too 
an exact value Ths = \/3/7r = 0.5513. 


~ 0.5405, ruling out the conjecture for 


• The conjectured universality of certain ghost coefficients seems to be somewhat prob¬ 
lematic in view of our data. While most coefficients come strikingly close to the conjectured 
values, the coefficient of c_i|0) deviates by 2% from the expected answer. Our predictions 
for the L = oo string field are likely to have a smaller error. 


• The ‘quasi-pattern’ observed in of an approximate factorization of LoTsiegei, where 
Tsiegei is the Siegel-gauge tachyon condensate, gets worse in our L = oo extrapolations. 
So this is definitely not an exact property. 


The correspondence between OSFT equations of motion in Siegel gauge and a certain 
equation for ghost number zero string fields [37| is still rather well obeyed, but does not 
improve in the L = oo extrapolations. This is likely an approximate pattern, as already 
suspected in [^|. 


We hope that our accurate data will stimulate new imaginative approaches to the problem 
of finding an exact solution. In the near future it will be possible to extract from our results new 
information about the kinetic term around the tachyon vacuum. It will also be straightforward 
to extend the methods of this paper to the computation of more general classical solutions of 
OSFT, which should provide more analytic clues. 

A more complete set of numerical data that it is practical to reproduce here (coefficients at 
higher levels, more significant figures etc.) will be be made available on-line p0| . 

To make this paper self-contained, we begin in the next section with a review of some basics. 


2 OSFT and the Universal Tachyon Condensate 


In this section we describe the basic setup for classical equations of motion in OSFT, with 

an emphasis on the symmetries obeyed by the tachyon condensate string field in Siegel gauge. 

While none of the results presented in this section are really new, we take the opportunity to 

spell out some facts - in particular our summary equation (2.13) - that are not widely known. 

action. By contrast, in t he ( L, 2L) scheme one keeps only the cubic terms in the action whose total level is < 2L. 
As discussed in section (4.3), (L, 3L) results display a much smoother dependence on L. 
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2.1 The Tachyon Condensate 

The action of OSFT takes the well-known (deceptively) simple form 

-S = -^ (^(^, Qb^) + ^ (T, T * . (2.1) 

This action describes the worldvolume dynamics of a D-brane specified by some Boundary CFT. 
The string field T belongs to the full matter+ghost state-space of this BCFT. In classical OSFT, 
T has ghost number one0. According to Sen’s conjecture |l|, the classical OSFT eom’s 

= 0 ( 2 . 2 ) 

must admit a Poincare invariant solution T = T corresponding to the condensation of the open 
string tachyon to the vacuum with no D-branes. The tachyon potential V(T) is given by @ 

V(l') 11 1 

where M is the brane mass. The normalized potential /(T) is expected to equal minus one at 
the tachyon vacuum, 

/(T) = -l. (2.4) 


2.1.1 Universality 

A basic property of the tachyon condensate string field T is universality 

(2.5) 


where 


n,. 


= Span{L™ 


ji 


T™ „ 6-fci ... c_i^ ... c_/JO) , ji >2,ki> 2, k >-1} 


( 2 . 6 ) 


with L™ denoting the matter Virasoro generators. The universal space is further decomposed 
into a direct sum of spaces with definite ghost number 


n 


univ 


®n£zT~i 


(n) 

univ ‘ 


(2.7) 


The restriction of the classical action to TCuniv can be evaluated using purely combinatorial 
algorithms that only involve the ghosts and the matter Virasoro algebra with c = 26 |Q, |^. It 
follows that the form of T does not depend on any of the details of the BCFT that defines the 
D-brane background before condensation. 

® Our conventions and notations are the same as Q. In particular we define the SL(2,R) vacuum |0) to have 
ghost number zero, and the ghost and antighost fields c{z) and b{z) to have ghost number one and minus one, 
respectively. 




2.1.2 Twist 


An obvious symmetry of the tachyon condensate is twist symmetry. The OSFT equations of 
motion admit a consistent truncation to twist even string fields |41], and indeed the tachyon 
condensate solution turns out to be twist even. In twist is defined simply as (—1)^°^^, 

so T contains only states with even level L = Lq + 1. 


2.1.3 Siegel gauge and SU(1,1) 

The Siegel gauge condition 6o'I' = 0 is particularly natural in level truncation since it is easily 
imposed level by level by simply omitting all Fock states containing the ghost zero mode cq. 
The Siegel gauge-fixed equations of motion 

Lo^ + bo{^ * ^) = 0 ( 2 . 8 ) 

admit a consistent truncation to the subspace of string fields which are singlets of SU(1,1) |^2|. 
The SU(1,1) symmetry in question is generated by 

oo oo oo ^ \ 

Q — 'y ) (^C—jibn b—nCn) X — ^ C—nCn) T — ^ ) f j (^•^) 

n=l n=l n=l ^ 

and the singlet subspace is defined as 


^ G Tisingl iff bo'if = g^ = X^ =Y^ = 0. 


( 2 . 10 ) 


Notice that acting on Siegel states Q is just ghost number shifted by one unit, so all states in 
'Hsingi have ghost number one. To show consistent truncation of equations (^) to the singlet 
subspace, we need to prove that if G TLsingh then 6o('I' =t='I') G 'Hsingh so that all components of 
'k outside TLsingi can be consistently set to zero. A simple argument is as follows. The generator 
A is a derivation of the *-algebra^, and commutes with 6o- Hence if A'h = 0, A 6o('k * 'k) = 0. 
Clearly Q is also zero on 6o('k * 'k), since ghost number adds under ^-product. By the structure 
of the finite-dimensionalQ representations of SU(1,1), a vector with zero Q and A eigenvalues 
must also have zero Y eigenvalue, that is, 6o('k * 'k) G TLsingh a-s desired. 

The SU(1,1) singlet subspace has a simple characterization in terms of the Virasoro genera¬ 
tors of the ‘twisted’ ghost conformal field theory of central charge —2 00, 


TLsingi = Span . . . L'^_fc„Ci|0) ,ki>2}0 TLmatter (2.11) 

®It is enough to notice that —2X = {Qs,co}, see ( |2.1(j ) below. Both Qb and co are derivations Q, and 
(anti)commutators of (graded) derivations are derivations. On the other hand, the generator Y is not a derivation. 

^Since SU(1,1) commutes with Lq, we can run the argument in the subspaces of Huniv with given Lo, which 
are finite-dimensional. 

® A proof of the equivalence of definitions 


(2.1C) and (2.11) for Hsingi can be found in [M, section 3. 
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L 

0 

2 

4 

6 

8 

10 

12 

14 

16 

18 

20 

Ml,i 

1 

2 

6 

17 

43 

102 

231 

496 

1027 

2060 

4010 

Nl 

1 

3 

9 

26 

69 

171 

402 

898 

1925 

3985 

7995 

Ml,2 

0 

1 

4 

12 

32 

79 

182 

399 

839 

1700 

3342 

N'l 

0 

1 

5 

17 

49 

128 

310 

709 

1548 

3248 

6590 


Table 3: Dimensions of some relevant subspaces of TCuniv 


where 

L'\ = Ll + kjf + 6k,0 

The statement that 


+ 00 

{k-n) : bnCk-n ■ ■ 

n=—oo 


( 2 . 12 ) 


Tc 


Siegel 


e K 


( 1 ) 


"Hsingl — 


Span{L'5_,^ ... L'3_k^ ci|0) , fc, > 2, j- > 2, Y^i + E-?* ^ 

summarizes all the known linear symmetries of the Siegel gauge tachyon condensate, 
exact constraints (quadratic identities |Q) are considered in section 


(2.13) 

2n} 

Other 


2.2 Level-Truncation and Gange Invariance 

We measure the level L of a Fock state with reference to the zero momentum tachyon ci|0), which 
we define to be level zero, in other terms L = Tq+I. As usual, the level truncation approximation 
(L, N) is obtained by truncating the string field to level L, and keeping interactions terms in 
the OSFT action up to total level N, with 2L < N < 3L. In our numerical work we have 
systematically explored both the (L, 2L) scheme, which is (naively) the most efficient, and the 
(L, 3L) scheme, which is the most natural. In section f4.3| we discuss some empirical differences 
between these two schemes. 

The most economic representation of Tsiegei is using the basis (|2.13|) , but unfortunately we 
have not found a simple algorithm to perform computations within the SU(1,1) singlet subspace^. 
We shall work instead with the universal basis ( |2.6|) using fermionic ghost oscillators. In this 
basis, the number of modes in Tsiegei truncated at level L (with L an even integer) is given 
by 

L/2 

iVL = E^2yi, (2.14) 

j=0 

where Mi,g denotes the number of Siegel Fock states in TCuniv with level I and ghost number g. 
®The twisted ghost Virasoro’s L'^ do not have simple conservation laws on the cubic vertex. 


10 























Mi^g which is computed by the generating function 


1 


I g P~‘^ ^ ^ ^ 


(2.15) 


The Siegel gauge-fixed eom’s (|2.8|) truncated at level L are a system of iV^ equations in iV^ 
unknowns. As we discuss in appendix A, the solution can be found very efficiently using the 
Newton method. By construction the resulting string field solves the truncated Siegel 


gauge eom’s with extremely good accuracy. However the full gauge invariant eom’s (^) impose 
an extra set of constraints on the solution. Recall that the BRST operator can be written as 


Qb = cqLo — 2boX + Q , 


(2.16) 


where 


Q= T. 


m — n 


Cm.O, 


ib—m—n ^ ^ 


m,n7^0 

m+riT^O 


n/0 


The extra conditions on a Siegel string field are then 


(2.17) 


QT 6oCo(^' * T) = 0. 

At level L, this equation entails extra constraints on 'T^g^g/, with 

L/2 

N’l = Y.M2j,2. 

1=0 

Table shows the numbers Ml,i, Nl, Ml ,2 and N'l up to L = 20. 


(2.18) 


(2.19) 


The role of equation ( 2.18 ) is simply to enforce extremality of the solution along gauge 
orbits. However, in principle there could be an issue about the non-perturbative validity of the 
Siegel gauge condition (are gauge orbits non-degenerate at the non-perturbative Siegel gauge 
vacuum? |^). Moreover, the level truncation procedure explicitly breaks gauge invariance, 
which is formally recovered only as L —> 00 . Thus equation (|2.18 ) gives an independent set of 
constraints which are not a priori satisfied by the level-truncated solution. If Siegel gauge is 
a consistent gauge choice and if gauge invariance is truly recovered in the infinite level limit, 
then we expect ( p.l8 ) to hold asymptotically as L —> 00 . This is a very non-trivial consistency 
requirement on TA 


Siegel ‘ 


Numerical evidence for this is examined in section 3.2. 


3 The Level-Truncated Tachyon Condensate 

Using the numerical methods outlined in appendix A, we have determined Tgiegei up to L = 18, 
both in the (L, 2L) and in the (L, 3L) schemes. {L, 3L) results appear to be better behaved (we 
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L = 4 

L = 6 

L = 8 

L = 10 

ci|0) 

0.548399 

0.547932 

0.547052 

0.546260 

c i|0) 

0.205673 

0.211815 

0.215025 

0.216982 

L™2Ci|0) 

0.056923 

0.057143 

0.057214 

0.057241 

C-3|0) 

-0.056210 

-0.057392 

-0.057969 

-0.058290 

6 2C 2Ci|0) 

-0.033107 

0.034063 

0.034626 

0.034982 

6 3C iCi|0) 

0.018737 

0.019131 

0.019323 

0.019430 

L™2C-i|0) 

-0.0068607 

-0.0074047 

-0.0076921 

0.0078698 

T™4Ci|0) 

-0.005121 

-0.005109 

-0.005102 

-0.005095 

L™2^-2Ci|0) 

-0.00058934 

-0.00062206 

-0.00063692 

-0.00064553 



L = 12 

L = 14 

L = 16 

L = 18 

ci|0) 

0.545608 

0.545075 

0.544637 

0.544272 

c i|0) 

0.218296 

0.219236 

-0.219942 

-0.220491 

i™2Cl|0) 

0.057252 

0.057256 

0.057257 

0.057257 

C-3|0) 

-0.058489 

-0.058625 

-0.058721 

-0.058794 

6 2C 2Ci|0) 

0.035225 

0.035402 

0.035535 

0.035640 

6 3C 1C1 0) 

0.019496 

0.019542 

0.019574 

0.019598 

L™2C-i|0) 

0.0079906 

0.0080782 

0.0081445 

0.0081966 

i™4Cl|0) 

-0.005090 

-0.005086 

-0.005082 

-0.005079 

T™2^-2Ci|0) 

-0.00065124 

-0.00065532 

-0.00065839 

-0.00066081 


Table 4: (L,3L) level-truncation results for the lowest modes of Tsiegei- 


come back back to this point in section O), and in this section we only consider this scheme. 
It is clearly impossible to reproduce here all the coefficients of the tachyon condensate up to 
level 18. We give some sample results in Table Q. Our complete numerical data will be made 
available at [|^]. 

In this section we perform some consistency checks of the level-truncation results, verifying 
some exact properties that the tachyon condensate must obey. 


3.1 SU(1,1) invariance 

We have systematically checked that our solutions for the tachyon condensate can be written 
in the basis ( 2.13| ), and thus obey the full SU(1,1) invariance. This property holds with perfect 
accuracy (that is, with the same precision as the number of significant digits that we keep, which 
is 15 for double-precision variables in C-I--I-). This is nice, but not surprising, since the SU(1,1) 
generators commute with Lq, and thus SU(1,1) is an exact symmetry of the level-truncated 
theory. 
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L = 6 

L = 14 

L = 00 

C 2Ci|0) 

0.00841347 

0.00257255 

-0.0000400232 

C- 4 C 1 IO) 

-0.0103276 

-0.00307849 

0.0000536768 

C iC 2|0) 

0.0107901 

0.00483115 

0.000005367 

L™ 2 C- 2 Ci| 0 ) 

0.000892329 

0.000612637 

0.00000877198 

L™3C iCi|0) 

-0.00212947 

-0.000877716 

0.00000163665 

C 6Ci|0) 

0.0130217 

0.00341282 

0.0000208782 

C 4C i|0) 

-0.0110576 

-0.00431119 

-0.000160066 

C- 3 C -2 0) 

0.00360400 

0.00160614 

-0.0000134344 

6 2C 3C iCi|0) 

-0.00306293 

-0.000919219 

-0.0000799493 

6 3C 2C iCi|0) 

-0.00324329 

-0.00114819 

-0.0000488214 

L™2'^-4'^1 0 ) 

0.000132483 

-0.000183042 

-0.0000162206 

L™2C-2C-i|0) 

-0.00188148 

-0.000811710 

-0.0000098375 

T™3C-3Ci|0) 

0.000834397 

0.000303847 

-0.0000004570 

L?4C 2Ci|0) 

0.000127107 

0.0000135260 

0.0000021124 

L™ 2 T™ 2 C- 2 Ci| 0 ) 

-0.000179524 

-0.0000980000 

-0.0000014704 

L™5C iCi|0) 

0.000903154 

0.000310410 

0.0000131051 

L™ 3 L™ 2 C-lCl| 0 ) 

0.000271962 

0.000105286 

0.0000014747 


Table 5: Sample (L, 2>L) level-truncation results for the out-of-Siegel equations of motion. The 
table shows data for L = 6 and L = 14, and L = oo extrapolations obtained from the data for 
2 < L < 14 with a polynomial fit in 1/L. 


3.2 Out-of-Siegel Equations 


We now turn to the crucial check of the extra conditions imposed by the full equations of motion 
before gauge-fixing. We were able to carry out this computation up to L = 14. Table ^ shows 
some sample results for the string field ( |2.18| ) evaluated for T = Tsiegei- The extra constraints 
are satisfied already very well at L = 6 , and significantly better at L = 14|3. This is happening 
thanks to large cancellations between the two terms in (2.18)'^, as can be easily checked by 
applying the operator Q to the results in Table 

Even more remarkable are the extrapolations of the data to L = oo, which give values two or 
three orders of magnitude smaller than the L = 14 results! Our extrapolation method consists 
in interpolating the data with a polynomial in 1/L of maximum degree (that is, with as many 
parameters as the number of data points). For example, for the mode c_ 4 Ci| 0 ) we have six 
data points (L = 4, 6 , 8,10,12,14) and we use a polynomial in 1/L of degree five. Empirically, 
this method gives better results (L = oo extrapolations closer to zero) than making fits with 
polynomials in 1/L of lower degree. 

This analysis leaves little doubt that the full equations of motion are satisfied as L —> oo. 


This behavior is common to the higher ievei modes not reproduced in Tabie ^ 

At L = 14, each term in (2.18) is typically one or two orders of magnitude bigger than their sum. 
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T m 
^2 

r m 

j m 

r m 

T m 
-^10 

T m 
^12 

J m 

^14 

T m 

2 

1.127927 








4 

1.069643 

1.079864 







6 

1.046467 

1.051898 

1.053517 






8 

1.034587 

1.037554 

1.040767 

1.036977 





10 

1.027439 

1.029304 

1.031367 

1.033082 

1.025628 




12 

1.022688 

1.023975 

1.025369 

1.026797 

1.027437 

1.017346 



14 

1.019312 

1.020257 

1.021261 

1.022317 

1.023271 

1.023102 

1.011026 


16 

1.016795 

1.017520 

1.018279 

1.019076 

1.019875 

1.020461 

1.019662 

1.006039 

oo 

0.999916 

0.999877 

1.00429 

1.00526 






Table 6: {L, 3L) level-truncation results for Schnabl’s quadratic matter identities. The table 

shows the values for the ratios Rn of equ. (^). 


3.3 Exact Quadratic Identities 


As pointed out by Schnabl |^], any solution of the OSFT eom’s must obey certain exact 
quadratic identities that follow from the existence of anomalous derivations of the star product. 
An infinite set of identities is obtained from the anomalous derivations — L'^ 2 n- They 

are 

^ 

= i-ir - ('hlcoLolT), (3.1) 

where T is a solution in Siegel gauge. 

In Table ^ we show the level-truncation results for the ratios 

54 (T|coL5;|T) 


Rn = (-1)^ 


65 (T|coLo|T) 


(3.2) 


which are of course predicted to be exactly one. The quadratic identities are satisfied quite well 
already at low levels, and the extrapolations to L = oo (performed again with polynomials in 
1/L of maximum degree) give really good results. 


Both the quadratic identities just analyzed and the out-of-gauge eom’s (2.18) are exact 
constraints on the solution that are broken by level-truncation. We have found that the level- 
truncated answers for this class of observables are very accurately converging to their known 
exact values as L —> oo. This is strong evidence for the idea that level-truncation is a convergent 
approximation scheme. We have also learnt that maximal polynomials in 1/L give very precise 
extrapolations. It seems safe to assume that this should be a universal feature, and in the 
following we shall adopt the same extrapolation technique to quantities whose exact asymptotic 
value is a priori unknown^. 

^■^Polynomials in 1/L have also been used in the extrapolation procedure of of [^. It was also noted in 
that large level results appear to have corrections of order 1/L, although there the definition of level is somewhat 
different. 
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4 Extrapolations to Higher Levels 


Encouraged by the successful extrapolations to L = oo described in the previous section, and 
inspired by Taylor’s analysis |^, we have set up a systematic scheme to extrapolate to higher 
levels the results for the vacuum energy and for the tachyon condensate string field. In this 
section we focus on the results for the vacuum energy, while in the next we shall examine the 
results for the individual coefficients of the tachyon condensate. 


Unless explicitly stated otherwise, the use of the (L, 3L) scheme is implied in the rest of the 
paper. We justify this choice in section 4.3, where we briefly contrast {L,2L) versus (L,3L) 
results. 


4.1 Extrapolations of the Tachyon Effective Action 


The basic idea of our method has already been explained in the introduction. The first step is 
the computation of the tachyon effective action Vl{T), obtained by integrating out the higher 
modes up to level L. Some details of how this is done numerically are explained in appendix A. 
Figure O shows the plots of Vl{T) for 0 < L < 16. There is good convergence as L increases, 
indeed the curves for L > 6 are indistinguishable on the scale of Figure 0. 


For our extrapolations, we focus on a interval for the tachyon vev around the non-perturbative 
vacuum. We take 0.54 < T < 0.55. The function V^\t), where M is an even integer < 16, is 
the extrapolation ‘of order M’ of the tachyon effective action at level L, and is constructed as 
follows. We fix the dependence on L by writing 




M/2+1 

E 


n=0 


an{T) 

{L + lY ’ 


(4.3) 


for some coefficients functions an{T). The functions an{T) are determined by imposing the 
conditions 


vI^\t)=Vl{T), forL = 0,2,---,M, V T E [0.54,0.55] . (4.4) 


In other terms, we interpolate the M/2 + 1 values {Vl{T)\L = 0, 2,... M} with a polynomial in 
1/(L + 1) that passes through all the data point^. 

Our best estimate for the tachyon effective action at level L is the function Figure 

1^ shows the plots of (T) for L between ten and infinity. The position of the minimum in 
each curve defines our order M = 16 estimates and for the tachyon vev and vacuum 
energy at level L. We can follow very clearly the behavior of the minima as L increases. The 

^^The rationale for using polynomials in 1/(1/ + 1) rather than 1/L is that we wish to include also the data 
for L = Q. This works somewhat better than excluding the L = 0 point and making extrapolations in 1/L. 
Committed readers can find more about this technicality in footnote H. 
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Potential 



Figure 2: Plots of the tachyon effective potential Vl{T) at level L, for L between zero (uppermost 
curve) and 16 (lowermost curve). The curves for L = 6,10,12,14,16 appear superimposed in 
the figure. 


energy falls below -1, reaches its lowest point in L = 28 curve, and then turns back to approach 
asymptotically the value = —1.00003 ! In Figure ^ we see the same phenomenon in a plot 
of as a function of 1/L. 

It is interesting to consider how the extrapolations change as we vary M. Table shows 
the estimates up to L = 18, while Table 12 (appendix B) shows the analogous estimates 

for the tachyon vev. By construction, the diagonal entries and are simply the exact 
values obtained by direct level-truncation at level L. One can observe from the tables that the 
method has remarkable predictive power. For example, by only knowing level-truncation results 
up to level 10, one can obtain the prediction = —1.0003780 for the energy at level 16, to 
be compared to the exact value Eiq = —1.0003755. We thus feel quite comfortable in trusting 
the extrapolations even for L large. Figure || and Table (already discussed in the introduction) 
illustrate the main features of the larger L results for the vacuum energy, for various M’s. It is 
pleasant to observe that, as M increases, the estimates E^^ (Table Q) quickly reach 

stable values, while E^'^ steadily approaches minus one|^. 

Finally we would like to comment on how results change if instead of using a polynomial extrapolation in 
1/(1/ + 1) we use an extrapolation in 1/L (excluding the L = 0 point), or alternatively we keep the L — 0 point 
and use a polynomial in 1/(L + a) for some other a > 0. One finds that for M = 16 the differences among all 
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Potential 



Figure 3: Plots of the order 16 estimates V^^\t) for the effective tachyon potential, for some 
sample values of L > 10. The minimum of each curve is indicated by a black dot, which by 
definition has coordinates {T^^\ ). The isolated uppermost plot corresponds to L = 10. 

To follow the curves from L = 10 to L = oo, focus on the position of the minima: as L increases, 
the dot moves from right to left (i.e., the tachyon vev decreases). As L —> oo, the curves crowd 
towards an asymptotic function with minimum at = (0.5405, -1.00003). 


It is remarkable that extrapolations to higher levels work so well. The data have a smooth 
and predictable dependence on L, which is very well captured by polynomials in 1/L. This 
property was not a priori obvious, and indeed it appears to be true only for (L, 2>L) data, as we 


shall see in section 4.3 


4.2 Comparison with Straightforward Extrapolations 

The method just described appears to work remarkably well. To which extent does the success of 
the method depend on the sophisticated idea of extrapolating the functional form of the tachyon 

these schemes are very minor, even for a wide range of reasonable values of a (say 0.1 < a < 3). For M < 10, 
including the L — 0 data (and using 1/{L + a)) works somewhat better than excluding it (and using 1/L). For 
example the prediction Elf obtained excluding L — 0 differs by the exact value by an error of 0.0003, which is 
30 times bigger than for the prediction obtained including L = 0. All of this scheme-dependence is expected to 
disappear for M large, and indeed is already irrelevant at M = 16. 
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Energy 



(1 

Figure 4: Plot of E\^ as a function of 1/L. The black dots represent the exact values up to 
L = 18 computed by direct level-truncation (Table 1, (L, 3L) scheme). To first approximation, 
the curve in the figure is roughly a parabola: since the energy overshoots -1 at 1/L = 1/14 ~ 
0.07, we have a visual understanding of the position of the minimum of the energy around 
1/L = (l/14)/2 = 1/28 ~ 0.036. 


effective action? We can answer this question by considering the more straightforward procedure 
of simply extrapolating the values of the vacuum energies El (as opposed to the full functions 
Vl{T)). We define the ‘straightforward’ order M estimate E^'^ at level L by considering the 
data {El\0 < L < M}, and interpolating them with a polynomial in 1/(L -|- 1) of maximum 
degree. (This is in complete analogy with (^)). The results for L < 18 in the (L,3L) scheme 
are presented in Table |^, while in Table ^ we give the extrapolations to L = oo. We see that 
for M < 10, the more sophisticated method gives much more accurate predictions (compare 
Table ^ and Table for example e[^ , and the exact value L^g®^). However for M > 10 
there is no significant difference between the two procedure We also compared the results 
for the individual coefficients of the tachyon string field obtained with the two procedures, and 
found a very similar pattern. 

A comparison with the results of [ p^ , which in our language correspond to M = 10, shows that the accuracy 
of the perturbative method of [sot seems comparable with the accuracy of the straightforward extrapolation. For 
M = 10 our non-perturbative method based on the tachyon effective action appears to be more accurate. 
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L = 4 

L = 6 

L = 8 

L = 10 

M = 4 

-0.98782176 

-0.99546179 

-0.99850722 

-1.0000023 

M = 6 


-0.99517712 

-0.99798495 

-0.99930406 

M = 8 



-0.99793018 

-0.99918359 

M = 10 




-0.99918246 


L = 12 

L = 14 

L = 16 

L = 18 

M = 4 

-1.0008372 

-1.0013461 

-1.0016765 

-1.0019017 

M = 6 

-1.0000079 

-1.0004169 

-1.0006692 

-1.0008317 

M = 8 

-0.99982545 

-1.0001796 

-1.0003843 

-1.0005057 

M = 10 

-0.99982266 

-1.0001750 

-1.0003780 

-1.0004979 

M = 12 

-0.99982226 

-1.0001739 

-1.0003759 

-1.0004947 

M = U 


-1.0001737 

-1.0003755 

-1.0004938 

M = 16 



-1.0003755 

-1.0004937 

M = 18 




-1.0004937 


Table 7: Estimates for the vacuum energy obtained from extrapolations of the effective 

tachyon potential, at various orders M and for L < 18. Data in the (L, 3L) scheme. By definition, 
the diagonal entries coincide with the exact computation from direct level-truncation 

at level (L,3L) (Table 1). 

We conclude that with the sophisticated procedure one can achieve remarkable accuracy even 
for small M, where a naive extrapolation would work quite poorly. However if one is willing to 
perform level-truncation up to level 12 or above, the simpler extrapolation procedure is equally 
effective. 

4.3 (L, 3L) versus (L, 2L) 

All the extrapolations described so far are for results in the (L, 3L) scheme. We have investigated 
the data in the (L, 2L) scheme and concluded that their behavior as a function of L is not nearly 
as smooth: as a consequence, extrapolations to higher levels are less reliable. A glance at Table 
1^ and Figure || is sufficient to illustrate our point. We are comparing the ‘straightforward’ 
extrapolations of the vacuum energy to L = oo, for various M’s, obtained with data in the 
{L, 3L) scheme, with the analogous quantities in the (L, 2L) scheme. While the (L, 3L) data 
have a really smooth dependence on M and converge nicely to -1, the (L, 2L) data have a much 
more irregular behavior. A similar pattern is observed for extrapolations at finite L: estimates 
of exact results at level L < 18 are not nearly as accurate in the (L, 2L) scheme as they 
are in the {L, 3L) scheme. An analogous behavior is found in comparing (L, 3L) and (L, 2L) 
data for the tachyon vev. We have also repeated for (L, 2L) data the full analysis based on 
extrapolations of the tachyon effective action, and found no improvement with respect to the 
straightforward extrapolations shown in Table p. 
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L = 4 

L = 6 

L = 8 

L = 10 

M = 4 

-0.98782176 

-0.99730348 

-1.0020443 

-1.0048888 

M = 6 


-0.99517712 

-0.99845611 

-1.0002959 

M = 8 



-0.99793018 

-0.99921882 

M = 10 




-0.99918246 


L = 12 

L = 14 

L = 16 

L = 18 

M = 4 

-1.0067851 

-1.0081397 

-1.0091556 

-1.0099457 

M = 6 

-1.0014693 

-1.0022814 

-1.0028762 

-1.0033304 

M = 8 

-0.99991100 

-1.0003187 

-1.0005753 

-1.0007448 

M = 10 

-0.99982332 

-1.0001767 

-1.0003811 

-1.0005023 

M = 12 

-0.99982226 

-1.0001738 

-1.0003758 

-1.0004946 

M = U 


-1.0001737 

-1.0003755 

-1.0004939 

M = 16 



-1.0003755 

-1.0004937 

M = 18 




-1.0004937 


Table 8: The estimates for the vacuum energy in the (L, 3L) scheme, obtained with the 

‘straightforward’ polynomial extrapolation in 1/(L + 1). 



(L,3L) 

(L,2L) 

II 

-0.998698 

-0.988625 

00 

II 

-0.999784 

-1.00261 

M = 10 

-1.00010 

-0.999316 

M = 12 

-1.00008 

-1.00048 

M = 14 

-1.00004 

-0.999655 

M = 16 

-1.00003 

-1.00057 


Table 9: The estimates for the asymptotic vacuum energy in the (L, 3L) and (L, 2L) 

schemes, obtained with the ‘straightforward’ polynomial extrapolation in 1/(T + 1). 

It would be interesting to explain these findings from an analytic point of view. The (L, 3L) 
scheme can be understood as a cut-off procedure in which only the kinetic term of the OSFT 
action is changed, such to give an infinite mass to modes with level higher than L. On the other 
hand, in the (L, 2L) scheme both the kinetic and the cubic term of the action are changed. This 
may explain why (L, 3L) data have a simpler dependence on L. 


5 Patterns of the Siegel Gauge Tachyon Condensate 


The individual coefficients of the tachyon condensate can also be extrapolated to higher levels. 
The more sophisticated method that we use is the following: We solve the classical equations 
of motion at level L and express all higher modes as functions of the tachyon vev T (see 
(A..12)). We then perform an extrapolation of these functions of T using a polynomial in 1/L 
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-0.9995 

- 



■. 





12 

14 

16 

-1.0005 

- 




'■ 

-1.001 

- 





-1.0015 

- 





-1.002 

- 





-1.0025| 

i 






Figure 5: Plot of the results in Table |9[ The continuous line represents the (L, 3-L) results, 
while the dashed line represents the {L, 2L) results. 


of maximum degree: this defines in the usual way the order M estimates T^^[r]. Finally by 


setting T = Tl we obtain extrapolations for the full tachyon string field. Tables and 
(appendix B) contain the extrapolations to T < 18 of the first two higher modes. In Table n 
we give the M = 16 results for the extrapolations to T = oo obtained using this method]^. 

We are finally in the position to look for analytic patterns in the tachyon condensate, in 


particular checking in a more reliable way the patterns observed in |^|. This was one of the 
motivations of our work. 


The Tachyon Vev 

Figure ^ is a plot of the estimate for the tachyon vev as a function of 1/L. The asymptotic 
value is ~ 0.5405. Clearly the conjecture ||2^ for an exact value Ths = is falsified]^. 

^®The straightforward method of simply extrapolating the coefficients with a maximal polynomial in 1/L is 
also possible: for M = 16 the difference with the data in Table ^ is very minor, at most or two units in the last 
significant digit. 

'^^This conclusion was already believed based on level-truncation results to level 10 but at the time of the 
Strings 2002 conference we had proposed that somehow this conjecture could be rescued. In an attempt to 
explain the puzzling overshooting of the vacuum energy in Table 1, we had suggested an ad hoc renormalization 
of the tachyon condensate by an overall multiplicative factor, such that the tachyon vev is exactly Ths- With 
this renormalization, the vacuum energy at L = 18 becomes very accurately -1. In view of the new results in 
and in this paper, clearly there is no need for any such mechanism. 
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Tachyon 



Figure 6: Plot of the estimate for tachyon vev as a function of 1/L. The dots represent 

the exact values for L = 2... 18. The dashed line represents the value Tjjs = \/3/7r. 


Universal Ghost Coefficients 


The most accurate pattern discovered in |37] is the remarkable universality of certain ghost 
coefficients. If one normalizes the tachyon mode to one, then the coefficients of the modes 
C-nb-mCi\0) (with n and m odd) appear to be the same for all known OSFT solutions. Assum¬ 
ing that this is the case, an analytic prediction for these coefficients is possible by looking at 
infinitesimal OSFT solutions corresponding to exactly marginal deformations of the BCFT. 


We reproduce our results for these coefficients for the tachyon condensate in Table 0 and in 
Table (this latter table is in appendix B). There is certainly a striking pattern. In particular, 
the results of Table 0 show that the pattern persists in the higher-level modes that were not 
explored in |^]. However the value for ri^i is puzzling. This result has not improved with respect 
to the L = 12 data already available in |37]0, and the 2% difference from the conjectured value 
would seem like a real onef^ 


In we performed an extrapolation of the data up to L = 12 using a fit a + b/L, and found rip = 0.411545. 
If one instead applies to the data up to L = 12 an extrapolation with a maximal polynomial in l/L (as advocated 
in this paper), one finds np = 0.415947. Changing extrapolation scheme or adding more levels does not seem to 
help in getting closer to the conjectured answer. 

^®Some of the higher-level coefficients (Table also have errors of a few percents, but this need not be a 
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L 

n,i 

r3,i 

^5,1 

r3,3 

rj,! 

r5,3 

2 

0.349483 






4 

0.375042 

-0.102499 





6 

0.386571 

-0.104743 

0.0547758 

0.0208371 



8 

0.393062 

-0.105966 

0.05544 

0.0209129 

0.0358037 

0.014195 

10 

0.397214 

-0.106707 

0.0558499 

0.0209927 

0.0361037 

0.0142012 

12 

0.400096 

-0.107201 

0.0561103 

0.0210499 

0.0363031 

0.0142291 

14 

0.402212 

-0.107553 

0.0562894 

0.0210917 

0.0364338 

0.0142521 

16 

0.403832 

-0.107818 

0.0564199 

0.0211232 

0.0365254 

0.0142701 

18 

0.405111 

-0.108023 

0.0565194 

0.0211478 

0.0365932 

0.0142843 

oo 

0.4160 

-0.1097 

0.05728 

0.02135 

0.03739 

.01456 

conj 

0.407407 

-0.109739 

0.0577148 

0.021328 

0.037483 

.01439 


Table 10: (L, 3L) numerical results for the pattern coefficients rn,m for the tachyon condensate 
string field, and their conjectured values. The L = oo results are obtained from the M = 16 
extrapolation procedure based on the effective tachyon potential. 


Although we hesitate to assign a precise error to our extrapolations to L = oo, it seems 
that this error should be smaller than than 2%. Indeed all the properties that we know for sure 
to be exact (the quadratic identities, the out-of-Siegel equations and of course the value of the 
vacuum energy), are obeyed with an accuracy of order 10“^ in the M = 16 extrapolations^. 


Lq T factorization 

In it was observed that the string field Lq Tsiegei is approximately factored into a matter 
times a ghost component. In view of our more precise data, we definitely conclude that this is 


just a rough pattern, as already believed in |^]. Indeed the pattern is seen to get worse in the 
L = oo extrapolations. Let us give a couple of examples. Assuming factorization, the normalized 
coefficient of L™ 2 C-i| 0 ) can be obtained by multiplying the two normalized coefficients of L™ 2 | 0 ) 


problem, since at higher levels we expect larger errors. 

^°While the energy is stationary and so it is affected quadratically by a small change in the coefficients, Schnabl’s 
identities and the out-of-Siegel equations vary linearly. It would take some conspiracy for the vev of of c_i|0) to 
have a 2% error, and at the same time the corresponding out-of-Siegel equation be so well obeyed (see the first 
line in Table which is linearly influenced by a small variation of this coefficient). 
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'^Siegel 

4 > 

L-2 

.1058 

0.1069 

L 4 

-.009343 

-0.009476 

L 2L 2 

-.001260 

-0.001221 

L-6 

.002648 

0.002691 

L-sL-s 

.0000135 

0.0000143 

1 

1 

.000575 

0.000594 

L- 2 L- 2 L -2 

-.000009 

-0.000016 


Table 11: The left column shows the normalized matter coefficients of the tachyon condensate 
string field for L = oo (from Table 13). The right columns shows the corresponding coefficients 
for the ghost number zero string field, obtained from an extrapolation to L = oo of data up to 
L = 14, with a maximal polynomial in 1/L. 


and c_i|0). In |^, using the numerical values at level (10, 30), this procedure gave a prediction 


with a 4% error; with the L = oo data in Table IJ, the error is increased to 8%. Similarly, 
applying the same procedure to L™4C_i|0), one finds only a 0.1% error at level (10,30), but a 


1% error using the values in Table 13 


Correspondence with Ghost Number Zero 


The last pattern noticed in [^| is a surprising coincidence between the OSFT equation and 
the equation 

(Lo - 1)^ + ^ * ^ = 0, (5.5) 

where <1? is a ghost number zero held. This equation can be solved in the universal space spanned 
by total Virasoro’s L_„ acting on the vacuum. It was found that the coefficients of the terms 
L-ki ■ ■ ■ L_fc^|0) in the solution of (|5.5|) are strikingly close to the normalized coefficients of the 


matter modes • • • L™^^ci|0) of Tsiegei- However the pattern is somewhat irregular and does 
not appear to systematically improve with level. This is conhrmed by our more precise data, 
a sample of which is presented in Table |^. There is no clear improvement with respect to the 


data in |^|, and we conhrm the conclusion that this is likely only a ‘quasi-pattern’. 


All in all, the patterns of |37| are still present in our more precise data, and one cannot help 
the feeling that they hint at some analytic clue. However we seem to conclude that all these 
patterns are not exact properties, with the possible exception of the universality conjecture for 
certain ghost coefficients. This last conjecture is generally well obeyed, but we found a disturbing 
discrepancy for the c_i|0) mode. 
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6 Concluding Remarks 


Our results support the idea that level-truncation is a completely reliable approximation scheme 
for OSFT, with a convergent limit as the level is sent to infinity. All available exact predictions 
(notably the value of the vacuum energy) are accurately confirmed by the data. No inconsisten¬ 
cies seem to arise from the fact that gauge-invariance is broken at finite level, indeed we found 
strong evidence that it is smoothly restored as L ^ oo. Quantities computed in level-truncation 
exhibit a predictable dependence on level which is very well approximated by polynomials in 
1/L (at least for the {L,3L) scheme). This allows reliable extrapolations to higher levels. Com¬ 
bining this observation with efficient computer algorithms based on conservation laws Q, we 
have developed very powerful numerical tools to study OSFT. 


In this paper we have focused on the universal subspace and obtained accurate data for 
the tachyon condensate. An obvious direction for further research is to use our results to learn 
about the kinetic term around the tachyon vacuum. The nature of this kinetic term is still 
rather mysterious. No perturbative open string states are expected to be present, and numerical 


evidence for this has already been obtained [31|. Our data will allow a more precise analysis, 
and hopefully give new analytic clues. 


The most intriguing aspects of the non-perturbative vacuum are related to the elusive closed 
string states. In OSFT, amplitudes for external closed strings (on a surface with a least one 
boundary) are given by correlation functions of certain gauge-invariant open string functionals 
||4^ , |4^ , [Till . It would be very interesting to compute such amplitudes in the non-perturbative 
vacuum. This should shed some light on the mechanism by which open string moduli are frozen 
in the tachyon vacuum, but closed string moduli are still present. There are promising ideas 
for how this may come about |H, ^], but the actual mechanism realized in OSFT is still 

unknown. 


Another avenue for future work is the application of our methods to more general classical 
solutions. It will be straightforward to extend our algorithms to include the matter states 
necessary to construct non-universal solutions, e.g. tachyon lump solutions and Wilson line 
solutions [^]. It would also be very nice to investigate numerically time-dependent solutions, 
and demonstrate the existence of tachyon matter [50| in OSFT. The study of several classical 
solutions will help to build the intuition that is needed for analytic progress. 

We hope that the precise data presented in this paper will encourage other physicists to think 
about analytic approaches to OSFT. Our analysis of the ‘patterns’ of the tachyon condensate 


throws some doubt even on the most robust conjecture proposed in |37], although it does not 
rule it out completely. The search for analytic clues in OSFT solutions is still very open. Yet 
we feel that our numbers for the tachyon string field must possess some hidden beauty, to be 
unveiled when an exact solution is found. 
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A The Numerical Algorithms 

In this appendix we explain some technical details about the algorithms that we have used to 
compute the universal star prodncts and find the tachyon solntion in level truncation. 

A.l Star Products from Conservation Laws 

The strategy for evaluating star prodncts using conservation laws is explained in detail in Q. 
Each Fock state in can be represented as a string of negatively moded ghost or Virasoro 

generators acting on the zero-momentum tachyon ci 10). The triple product of three such states 
is evaluated recursively by converting a negatively moded generator on one state space to a sum 
of positively moded generators acting on all three state spaces, 

(A_fc$i,T2,4'3) = (A.l) 

rfc( 4 >i, Ta, ^3) + Ta, ^3) + ($1, ^ Pn^n^ 2 , ^3) + (^1, ^2, ^ , 

n>0 n>0 n>0 

where A_k is any constructor symbol and the coefficients appearing in this ‘conservation law’ 
are computed from the geometry of the Witten vertex Q . All triple products in are thns 

reduced to the coupling (ci,ci,ci) of three tachyons. 

Once the triple prodncts are known, star products are easily obtained by inverting the non¬ 
degenerate bpz inner product. If {4'^} is a Fock basis for we define the dual basis {'h*} 

( 2 ) 

of by the bpz pairing 

(XI/,, 4 /^) =<5^-. (A.2) 

Then 

(A.3) 

k 

We automated this algorithm on a C-|--|- computer code. We briefly highlight some features 
of our implementation; 

• We use the factorization of the star prodnct into matter and ghost sectors. The algorithm 
to find the triple products is executed separately in the two subsectors. 
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• We use cyclic and twist symmetry of the vertex to reduce the computation to triple prod¬ 
ucts ('hi, 'I'j, 'I'fc) with a canonical ordering i < j < k. 

• While in the matter sector the algorithm can be implemented in a straightforward way, 
in the ghost sector we need to face a slight complication. We are ultimately interested in 
triple products of ghost number one states, but the use of fermionic ghost conservation 
laws necessarily brings us outside the ghost number one subspace. We found it most 
efficient to use only conservation laws for the b-k oscillators. A single application of a 
6-ghost conservation law reduces the evaluation of a (1,1,1) product (ghost number one in 
all three slots) to a sum of terms of the form (1,1,1) and (0,1,2). Products of this latter 
type can be treated by applying a 6-conservation law to the first state (of ghost number 
zero), obtaining a sum of terms (1,1,1) and again (after cyclic rearrangement) (0,1, 2). It 
is easy to show that this algorithm always terminates on the product of three tachyons. 
So we see that we only need to consider triple products of the form (0,1,2) besides the 
standard products (1,1,1). 


After each application of a conservation law, the resulting triple products on the r.h.s of 


(A.l) are processed using the Virasoro algebra or the ghost commutation relations, till 
all states are reduced to the Fock basis ( |2.6D with the canonical ordering ji > ... > jp, 
ki >...> kg, li >...> Ir- The evaluation of expressions like • • • L_„.ci|0), 

with ni > re 2 • • • > n* (and similarly for the ghosts) is thus a basic elementary operation. 
There is a critical gain in efficiency in evaluating beforehand all such expressions (up to 
the desired maximum level) and reading the results from a file, rather then re-computing 
them each time. The size of such a file grows only linearly with the number of modes, 
whereas the table of triple products grows cubically, so this strategy is not problematic 
from the point of view of memory occupation. 


This algorithm can be easily extended to evaluate more general star products of string 
fields belonging to a larger space than Ttuniv, for example the space relevant for tachyon lump 
solutions [^] or Wilson line solutions [25|. One needs to enlarge the algebra of matter operators 
and consider the appropriate conservation laws. 


A.2 Solving the Equations of Motion 

Once all triple products at level L have been computed, the evaluation of the star product of two 
Siegel gauge string fields at level L involves 0{N^) algebraic operations. It is clearly desirable to 
have an algorithm to solve the classical eom’s that requires as few star products as possible. We 
tried various options, which can all be represented as a recursive procedure = E(T^”)), 

where T = F’('I') implies that T is a solution. 
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The most obvious idea is to invert the kinetic term in Siegel gauge and define 

F(T) =(A.4) 
To 

Clearly this iteration cannot converge since T(AT) = A^-F('k). There is a simple way to fix this 
problem, defining 

where T[<1>] indicates the coefficient of ci|0) in the string field d>. Unfortunately the algorithm 
based on the recursion F still fails to converge, and generically falls into stable two-cycles. An 
improved recursion is 

Fo(T) = aT + (a - 1)—T * T (A.6) 

To 

where a is a real number which is chosen randomly in some reasonable interval (say 0.2 < a < 
0.8) at each iteration step. This randomization breaks the cycles and the algorithm converges 
to a unique solution in about 20-30 steps. (The algorithm stops when the eom’s are satisfied 
with the same accuracy as the accuracy of double-precision variables in which have 15 

significant digits). This algorithm is very robust with respect to the choice of the starting point 
To, in fact at any given level L we found only one non-trivial solution. 

A more efficient approach is the standard Newton algorithm. Recall that given a system of 
N algebraic equations in N variables, fi[xQ.] = 0, 1 < i,a < N, the Newton recursion is 

^in+l) ^ ^(n) _ (A.7) 

where the matrix Mia[x] is defined as 

= (A.8) 

OXa 

In onr case, the truncated Siegel equations of motion are a system of Nl algebraic equations 

in Nl variables (Table ^) and this method can be directly applied. It is interesting to write 

the Newton algorithm as a recursion for the Siegel string field itself. One finds the compact 
expression 

Here the operator Qijf is defined by 

T = QsTT * TT * T (A.10) 

for any ghost number one string field T. The inverse operator is naturally defined by 
projecting onto the Siegel subspace. In other terms, for any ghost number two string field S, 
we look for the ghost number one string field that obeys 

^0 (Q<t^5]) = 0, coboF, = coboQii,{Q^^)'E,. (A.11) 
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The operator Q-qt has a natural physical interpretation: If is a solution of the OSFT eom’s, 
then Qq, is the new BRST operator obtained expanding the OSFT action around 'h. Thus 
as we approach the fixed point of the Newton recursion, Qq,{n) becomes a better and better 
approximation to the BRST operator around the tachyon vacuum. 


In level-truncation, the action of the operator Qq, in the Siegel subspace is represented by 
an Ni X Nl matrix. Since there is an order 0{Nf) algorithm to invert a matrix, the Newton 
recursion is not significantly more time-expensive than the evaluation of a single star product. 
The Newton algorithm is very fast, effectively squaring the accuracy at each step, and the 
solution is reached in four or five iterations. On our pc, the complete algorithm (computing 
the vertices from scratch and hnding the tachyon solution) takes less than 10 seconds at level 
(10,20), and less than a minute at level (12,36)! There is however a rather critical dependence 
on the initial conditions: one finds convergence only from a starting point sufficiently close to 
the solution (it is enough to take e.g. = 0.5 Clio)). 


We compared the solution obtained with the Newton method with the solution found with 
the alternative algorithm described above, finding exact agreement up to level 16. (At level 18 
the recursion ( A.6| ) runs too slowly on our pc). This gives a strong check on the correctness of 
the solutions. As another check, we compared our results at level (10,20) with the results of 
Moeller and Taylor @0, finding agreement to the tenth significant digit. 


A.3 Tachyon Effective Action 

To compute the tachyon effective action, we write the string field as 


Tl = Tci|0) + Tl, 


(A.12) 


where contains all the modes up to level L, except ci|0). For a given numerical value of of 
the variable T, we solve the classical OSFT equations of motion for all the higher modes, using 
the Newton method. This gives ^l\T'] = Tci|0) -|- ^l[T] as a function of T. Plugging ^l[T] 
into the OSFT actionp^ we obtain the effective tachyon potential Vl{T). 

The Newton algorithm that finds the solution ^l[T] fails to converge if the variable T is 
outside an interval yve find for example ~ —0.1 and r{g“* ~ 0.7 (notice 

that both the tachyon vacuum and the perturbative vacuum are safely inside the convergence 
region). The failure of the numerical algorithm can be explicitly traced to the existence of other 


branches in the tachyon effective action. This phenomenon has been studied in |2C, 35], where 

it has also been related to the non-perturbative failure of the Siegel gauge condition. In this 

paper we only need Vl{T) in an interval around the non-perturbative vacuum, which we take 

thank the authors of for making their full results available to us. 

More precisely, Vl{T) = /(^_l[T]), where is defined in (2.3). 
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to be 0.54 < T < 0.55. We postpone a more detailed investigation of the global behavior of the 
tachyon potential. 

B Some Further Numerical Data 



L = 4 

L = 6 

L = 8 

L = 10 

M = 4 

0.54839904 

0.54849677 

0.54814406 

0.54777626 

II 


0.54793242 

0.54711284 

0.54639593 

M = 8 



0.54705245 

0.54626520 

M = 10 




0.54626093 


L = 12 

L = 14 

L = 16 

L = 18 

M = 4 

0.54745869 

0.54719393 

0.54697362 

0.54678893 

M = 6 

0.54581507 

0.54534684 

0.54496539 

0.54465027 

M = 8 

0.54561932 

0.54509452 

0.54466463 

0.54430807 

M = 10 

0.54560864 

0.54507682 

0.54464004 

0.54427703 

M = 12 

0.54560809 

0.54507524 

0.54463714 

0.54427267 

M = U 


0.54507515 

0.54463683 

0.54427204 

M = 16 



0.54463682 

0.54427198 

M = 18 




0.54427196 


Table 12: Estimates for the tachyon vev obtained from extrapolations of the effective 

tachyon potential, at various orders M and for L < 18. Data in the (L, 3L) scheme. By 
definition, the diagonal entries coincide with the exact computation from direct level- 

truncation at level (L,3L) (Table Q). 
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Table 13: Asymptotic values for the first coefficients of the tachyon condensate string field, 
compared with the L = 18 data. The L = oo results are obtained from the M = 16 extrapolation 
procedure based on the effective tachyon potential. Data in the {L, 3L) scheme. 
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L = 4 

L = 6 

L = 8 

L = 10 

M = 4 

-0.20567285 

-0.21119493 

-0.21392087 

-0.21552208 

M = 6 


-0.21181486 

-0.21499106 

-0.21690691 

00 

II 



-0.21502535 

-0.21697620 

M = 10 




-0.21698254 


L = 12 

L = 14 

L = 16 

L = 18 

M = 4 

-0.21656852 

-0.21730328 

-0.21784642 

-0.21826369 

II 

-0.21818110 

-0.21908711 

-0.21976325 

-0.22028663 

M = 8 

-0.21827982 

-0.21920964 

-0.21990503 

-0.22044413 

M = 10 

-0.21829559 

-0.21923573 

-0.21994128 

-0.22048993 

M = 12 

-0.21829570 

-0.21923600 

-0.21994171 

-0.22049051 

M = 14 


-0.21923603 

-0.21994180 

-0.22049069 

M = 16 



-0.21994181 

-0.22049069 

M = 18 




-0.22049069 


Table 14: Estimates for the vev of c_i|0), obtained from extrapolations of the effective tachyon 
potential, at various orders M and for L < 18. Data in the {L, 3L) scheme. By definition, 
the diagonal entries coincide with the exact computation from direct level-truncation at level 


(L,3L) (Table I). 



L = 4 

L = 6 

L = 8 

L = 10 

II 

0.056923526 

0.057062755 

0.057068423 

0.057045668 

M = 6 


0.057143493 

0.057209039 

0.057229308 

00 

II 



0.057214101 

0.057239805 

M = 10 




0.057241066 


L = 12 

L = 14 

L = 16 

L = 18 

M = 4 

0.057018045 

0.056991677 

0.056968053 

0.056947298 

M = 6 

0.057233609 

0.057231744 

0.057227479 

0.057222387 

M = 8 

0.057248895 

0.057251070 

0.057250189 

0.057247946 

M = 10 

0.057252093 

0.057256442 

0.057257742 

0.057257584 

M = 12 

0.057252005 

0.057256182 

0.057257253 

0.057256834 

M = 14 


0.057256190 

0.057257279 

0.057256887 

M = 16 



0.057257279 

0.057256886 

M = 18 




0.057256885 


Table 15: Estimates for the vev of L™2Ci|0), obtained from extrapolations of the effective 
tachyon potential, at various ‘orders’ M and for L < 18. Data in the (L, 3L) scheme. By 
definition, the diagonal entries coincide with the exact computation from direct level-truncation 
at level (L,3L) (Table 
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L = 10 

L = 12 

L = 14 

L = 16 

L = 18 

conj 


0.0259407 

0.0261085 

0.0262255 

0.0263037 

0.0263594 

0.027063 

r7,3 

0.0104886 

0.0104825 

0.0104948 

0.0105063 

0.0105159 

0.0105868 

r5,5 

0.00658638 

0.0065773 

0.00658224 

0.00658778 

0.00659265 

0.0066192 

Gi,i 


-0.0200117 

-0.0201181 

-0.0201948 

-0.0202469 

-0.0208326 

^9,3 


-0.00818159 

-0.00817378 

-0.00818 

-0.00818657 

-0.00824041 

r7,5 


-0.00525767 

-0.00524794 

-0.00524929 

-0.00525185 

-0.00526601 

ns,! 



0.0161045 

0.0161778 

0.0162318 

0.0167396 

ni,3 



0.00662999 

0.00662276 

0.0066262 

0.0066689 

?’9,5 



0.00431978 

0.00431134 

0.00431141 

0.00431915 

n,7 



0.00315946 

0.00315272 

0.00315241 

0.0031549 

ri5,i 




-0.0133607 

-0.0134141 

-0.013871 

ns,3 




-0.0055259 

-0.00551968 

-0.005553 

ni,5 




-0.00363327 

-0.00362626 

-0.003629 


Table 16: (L, 3L) numerical results for the pattern coefficients for the tachyon condensate. 
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